Relaxation Spectra at the Glass Transition: Origin of Power Laws 
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We propose a simple dynamical model of the glass tran- 
sition based on the results from a non-randomly frustrated 
spin model which is known to form a glassy state below a 
characteristic quench temperature. The model is character- 
ized by a multi-valleyed free-energy surface which is modu- 
lated by an overall curvature. The transition associated with 
the vanishing of this overall curvature is reminiscent of the 
glass transition. In particular, the frequency-dependent re- 
sponse evolves from a Debye relaxation peak to a function 
whose high-frequency behavior is characterized by a non- 
trivial power law. We present both an analytical form for 
the response function and numerical results from Langevin 
simulations. 
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In supercooled liquids, the glass transition is heralded 
by anomalously slow relaxations The phase below 
the glass transition temperature is non-ergodic and the 
dynamics is characterized by "aging" ^ . Recent experi- 
ments indicate that the the approach to the glass transi- 
tion has some universal features |^ when viewed in terms 
of the frequency-dependent response of the system. Some 
of these features are also shared by spin glass transitions 
1^. Theoretical research has focussed on spin-glass-like 
phenomena in nonrandom systems Q| and some of these 
models have been shown to be equivalent to mean-field 
spin-glass models |^,^. The dynamics of certain frus- 
trated XY models have also been shown to exhibit glassy 
behavior |^]. 

In this paper, we present a dynamical model of the 
glass transition and analyze the trends in its frequency 
dependent response. The model is characterized by a 
multi-valleyed free-energy surface modulated by an over- 
all curvature and the glass transition is associated with 
the vanishing of this overall curvature. A multi-valleyed 
free-energy surface has long been associated with glasses 
iQ. The features of the present model were deduced 
from simulations of a nonrandomly frustrated spin model 
1^. The frequency-dependent response of the system 
can be calculated in closed form, and its low tempera- 
ture approximation becomes exact in the high-frequency 
limit. The results show that, as the overall curvature 
approaches zero, the frequency-dependent susceptibility 
changes from a characteristic Debye relaxation peak to a 
power-law spectrum reflecting the distribution of curva- 
tures of the valleys. The predicted changes agree qualita- 



tively with the generic features identified in experiments 

We begin with a description of the results of simula- 
tions of the frustrated spin system. The model consid- 
ered is an Ising antiferromagnet on a deformable trian- 
gular lattice. Detailed analysis of this model has shown 
that there is a first-order transition from the disordered 
paramagnetic state to an ordered "striped" phase [§,0 . 
Frustration plays a crucial role in this phase transition 
|]l0| . The model is characterized by two order parame- 
ters: a staggered magnetization with Ising symmetry and 
a shear distortion with three-state Potts symmetry. 

Monte Carlo simulations have demonstrated that 
the time evolution of various physical quantities, such 
as the average energy, undergoes a qualitative change as 
the system is quenched below a characteristic tempera- 
ture T*, below the ordering transition To investigate 
the nature of the transition at T* , we computed the fluc- 
tuation metric, fl{t) ||ll|. This function has been shown 
to be sensitive to ergodicity breaking and has been used 
to study the glassy phase in Lennard- Jones systems . 
The metric is defined by: 

n{t) =< {iy{t) -vf> ; u{t) = - f dt'E{t') (1) 



Here i> is the ensemble average of the energy. If a system 
is ergodic, the fluctuation metric decays to zero at times 
longer than the equilibration time and behaves as 1/t 
for exponential relaxations \ ^a\ . Breaking of ergodicity is 
signaled by the appearance of a non-zero long-time limit 
and a trajectory dependence lllfl] . Fluctuation metrics 
obtained from our simulations (Fig. (|^)) show that the 
supercooled state at temperatures above T* is ergodic 
but the state below T* exhibits broken ergodicity. 

A free-energy function was obtained from the prob- 
ability distribution of energy Q assuming that quasi- 
equilibrium is reached in each of the valleys The 
resulting function is a one-dimensional projection of the 
actual, multidimensional free-energy surface. As seen in 
Fig. (|^c), at temperatures above T*, there is a well de- 
fined valley centered at an energy corresponding to zero 
shear ||] . The curvature of this valley decreases and sub- 
valley structures become more pronounced as T* is ap- 
proached and, at the ergodicity-breaking transition, one 
is left with a multi-valleyed free energy landscape. These 
observations are consistent with the idea of a shear in- 
stability, which is known to exist in this model pO|, oc- 
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curring at T*. The temperature T* separates a high- 
temperature regime where the valleys with finite shear 
distortions are not accessible from a low-temperature 
regime where these valleys become accessible. The loss 
of ergodicity is associated with divergence of the aver- 
age trapping time in these valleys and is reminiscent of 
weak-ergodicity breaking plp^] . 



tlme(]VICS) 
50000 



time (MCS) 

50000 lOOOOO 




Energy (E) 

FIG. 1. Monte Carlo results for r2"^(t);(a)r > T* and 
(b)T < T* . Fig. 1(c) depicts the structure of the free-energy 
surface for T > T* (solid line) and T < T* (dashed line). 
These results are from a system of size 96x96 

Based on these observations of the non-randomly frus- 
trated spin model, we have constructed a simple dynam- 
ical model for the glass transition. The dynamics is de- 
scribed by a Langevin equation based on a multi-valleyed 
Hamiltonian H{(j)) where (j) is an order parameter, similar 
to the shear distortion in our spin model. A particularly 
simple form which has the characteristics indicated in our 
simulations (c/Fig. (|l]c)), consists of piecewise parabolic 
valleys within a megavalley characterized by an overall 
curvature. The distribution of curvatures and depths of 
the subvalleys, together with the overall curvature char- 
acterize the model completely. The glass-transition is 
associated with the vanishing of the overall curvature. 

The Hamiltonian describing the piecewise parabolic 
model is: 

= £ H„^i„, (2) 



where {fJ-n} is the set of characteristic functions on the 
non-overlapping segments [0° — A„; -f- A„], covering 
the domain of (t>(t), and c„ are chosen such that H[(j)\ 
is continuous. The parameter R determines the overall 
curvature and the individual well curvatures are deter- 
mined by r„. The relaxational dynamics of the model is 
described by the Langevin equation 



with the Gaussian noise 

< vit)v{t') >=TS{t-t'). 



(3) 



(4) 



We are interested in studying the frequency dependent 
response of this system as the overall curvature is tuned 
to zero. The response can be calculated by constructing 
an appropriate generating function. Integrating over the 
Gaussian noise, leads to the dynamical action |13] 



dt 



d 



5H[cj>] 



dV" 5(t>{t) 
which defines the generating function 



(5) 



z[h] = j vm 
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dt 



exp|-^5[<^]+ J 



dth{t)(l)(t) 



(6) 



For an ergodic system, the response can be related to a 
correlation function through the fluctuation dissipation 
theorem. Approaching the glass transition temperature 
from above, we can use this relationship to calculate the 
imaginary part of the frequency-dependent susceptibility 
X{lu) from a knowledge of the correlation function x 

x^^<m> ■,<m>^^^^z[h]- 
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Introducing the "centered" order parameter, Xn — 4> ~~ 
(j)^ and "effective" subwell curvature Rn = rn + R, and 
using the non-overlapping property: ^mlJ'n = Smn/J^n, the 
dynamic action can be rewritten as: 



^ dtY, 



n=0 



/ d 

(^— + i?„)2:„(i) 



(7) 



The non-overlapping property also leads to the generat- 
ing function being written as a sum of functions belong- 
ing to each segment: Z[h\ = X]n=o°' Zn[h], where the 
subwell generating functions are given by: 



Zn[h] 



P.„.„ll.„||exp|-^ 



dt( 




In the above equation, we have introduced the variables: 
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If the subwells were not constrained to exist in finite re- 
gions of the order parameter space, then would 
be a constant and the response of the system obtained 
from the resulting generating function would be a sum 
of Lorentzians. It will be shown below that this is very 
nearly the case in the limit of large frequency (scale set 
by the noise strength F). 

Because of the non-overlapping property (c/Eq. (||)), 
the susceptibility can be expressed in terms x„ and 0„ 
which denote the susceptibility and average order param- 
eter of the individual, decoupled wells. The calculation 
of Xn defies any simple analytic approach, because of 
the finite region of functional integration in the partition 
function. However, if the noise is small enough and the 
system never escapes its subwell, we can neglect all these 
effects of finiteness of the wells. It can be shown that 
the error vanishes as we increase ^„ = i^n^n/^r. In the 
zero noise limit all [h] are the same and do not depend 
on h. For finite noise, the same effect can be achieved 
by increasing the frequency associated with the response 
function and therefore, the low-noise approximation is 
less restrictive for high frequencies. 

The simplest case is when the overall curvature, R van- 
ishes. In this limit, which we associate with the glass 
transition, the frequency-dependent response Xn is given 
by: 



(9) 



Here, uj^ = m'^+Rn^. The "effective" susceptibility x =< 
Xn > is given by the weighted sum over all the subwells: 



Pn 



'z' 



Since Xn depends on n only through the curvature r„ 
{R — 0), the susceptibility x can be expressed in terms 
of the probability of encountering a given curvature: 



X 



dr- 



P{r) 



.Pir) 



PnS{r - r„) (10) 



In general, P(r) depends on the distribution of curva- 
tures and the 'time' spent in a well with given curvature. 
The latter is determined by the depth of the well and 
the overall curvature. However, in the low-noise regime, 
where we have argued that p„ is independent of n, P{r) 
reflects the 'intrinsic' distribution of curvatures in the 
model Hamiltonian. If this distribution is described by 
a power law: r", then within our approximation which 



is exact in the high-frequency limit, the response of the 
system reflects this power law. 



dr- 



OJ 



(11) 



When the curvature R is non-zero, the probabilities, 
Pn, get multiplied by exp(— i?^(/)°^/2F) and even for ex- 
tremely small overall curvature, R, the distribution P{r) 
changes considerably due to a different sampling of the 
wells. As R increases only subwells centered close to 
= are sampled appreciably and, in the limit of large 
i?, the response of the system is dominated by a single 
parabolic well (Debye response). The intrinsic distribu- 
tion of curvatures is reflected in the susceptibility only in 
the limit of i? = 0. This result shows that by tuning the 
overall curvature of a multi-valleyed free-energy surface 
the high-frequency response of the system can be changed 
from a pure Debye spectrum (x — w~^) to a nontrivial 
power law which reflects the distribution of curvatures of 
the subwells. 
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FIG. 2. Plots of X{lj). (a) R = 0.001, increasing values of 
F; O.Ol(solid), O.l(dotted), 0.5(dashed); (b) F = 0.01 increas- 
ing values of R; O(solid), 0.001 (dashed), 0.01(long-dashed);(c) 
F = 0.5, at different values of R; O.OOl(dot-dashed), 
0.05(solid), O.l(long-dashed), 0.2(dashed); and (d) Suscepti- 
bility at i? = and F — 0.01 for two different values of a. 
The straight line has a slope of —0.5. 

To demonstrate these features explicitly and to analyze 
the effect of large noise, we performed simulations of the 
Langevin equation (|3|). The subwells were chosen to be 
distributed uniformly (a = 0) and the largest subwell 
curvature was determined by considering the numerical 
stability of the system which required that in each well 
(j)At < 2A„ In each run the starting point was near the 
bottom of a randomly chosen subwell. We calculated the 
susceptibility as an average over a large number ('~ 1000) 
of runs and each run was ~ 10000 Monte Carlo Steps 
long. The results are shown in the Fig (||). 

For zero overall curvature, irrespective of the noise 
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strength, the susceptibihty reflects the distribution of the 
subwefl curvatures (Fig. (||)d). 

For a non-zero overall curvature, one can distinguish 
between two distinct regimes (a) low noise such that there 
is no escape from wells and (b) noise large enough to 
make escape possible The behavior in the (a) regime 
is summarized in Figs. (|^a) and (^). At these noise 
strengths, a finite value of the overall curvature leads to 
a large damping of the contributions from all subwells 
centered away from ^ = 0, and one observes a superpo- 
sition of very few Lorentzian peaks (Fig. (|2fD)). If the 
noise strength is increased, at small i?, the susceptibility 
starts to pick up contributions from all the subwells and 
the response becomes non- Debye like (c/Fig. (|a)). 

The response of the system is qualitatively different in 
regime (b), as illustrated in Fig (^). The finiteness of 
the wells becomes a crucial factor and, at some minimum 
value of R, escape from the most distant wells became 
possible. The dominant contribution to the susceptibil- 
ity arises from the relaxation towards </) = and at large 
R there is a crossover to a single Debye peak with no 
evidence of the internal structure. In the intermediate 
regime, the susceptibility reflects a superposition of both 
probing the distribution and relaxation towards the cen- 
ter. In this case the shape of the curve is sensitive to the 
details of all the distributions: centers of subwells, their 
widths (depths) and curvatures. 

The numerical and analytical results demonstrate that 
by tuning R or the noise strength F, one can reach a 
special "glassy" phase with a frequency-dependent sus- 
ceptibility obeying a non-trivial power law (Eq. ( pi] ) . For 
large noise and curvature, one observes a Debye peak de- 
termined by the overall curvature R. Starting from this 
liquid-like phase, the power-law regime is reached by de- 
creasing the curvature. In the low-noise limit (regime 
(a)), the phase with the Debye-like response resembles 
a frozen amorphous or polycrystalline state since there 
is no escape from subwells. The location of the Debye 
peak is determined not by the overall curvature R but 
by the average curvature of the subwells located near 
(j) = 0. This frozen phase is unusual since the system 
still relaxes in a multitude of subwells and is far from 
equilibrium. The frozen phase can be made to approach 
the "glassy" phase by increasing the noise strength. The 
power-law behavior characterizes the glass and is not ob- 
served in either the liquid or the frozen phase. 

The transition from the "liquid" to the "glassy" phase, 
as depicted in Fig. (^s) is similar to experimental obser- 
vations in supercooled liquids dJ^. In both theory and 
experiment, the shift of the Debye peak towards lower 
frequencies is accompanied by a change in the power law 
which characterizes the high-frequency response. The 
scaling features observed in experiments |^ are intriguing 
and probably reflect some interesting physical features of 
glass formers. This physics would dictate the type of 
distributions of curvatures and heights which enter our 



theoretical model. As mentioned before, all these distri- 
butions are crucial in determining the behavior of x{^) 
in the regime intervening between the large R and R — 
phases. We have made no attempt to tailor these dis- 
tributions to reproduce the scaling behavior observed in 
experiments. 

In conclusion, we have shown that our model of the 
glass transition predicts changes in the high-frequency 
response of a system approaching the glass transition. 
This model can be taken only as a schematic picture of 
a real glass transition but suggests that a multi-valleyed 
free-energy surface with an overall curvature can be taken 
as a useful template for glassy systems. The transition 
associated with the vanishing of the curvature is a new 
type of critical point and our analysis shows that this 
"phase transition" is characterized by a non-trivial, high- 
frequency response. 
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